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ABSTRACT 

We provide a detailed theoretical study aimed at the observational finding about the 
v Octantis binary system that indicates the possible existence of a Jupiter-type planet 
in this system. If a prograde planetary orbit is assumed, it has earlier been argued 
that the planet, if existing, should be located outside the zone of orbital stability. 
However, a previous study by Eberle & Cuntz (2010) [ApJ 721, L168] concludes that 
the planet is most likely stable if assumed to be in a retrograde orbit with respect to 
the secondary system component. In the present work, we significantly augment this 
study by taking into account the observationally deduced uncertainty ranges of the 
orbital parameters for the stellar components and the suggested planet. Furthermore, 
our study employs additional mathematical methods, which include monitoring the 
Jacobi constant, the zero velocity function, and the maximum Lyapunov exponent. We 
again find that the suggested planet is indeed possible if assumed to be in a retrograde 
orbit, but it is virtually impossible if assumed in a prograde orbit. Its existence is found 
to be consistent with the deduced system parameters of the binary components and of 
the suggested planet, including the associated uncertainty bars given by observations. 

Key words: methods: numerical — methods: 3-body simulations — techniques: 
Maximum Lyapunov Exponent — stars: individual (y Octantis) 



1 INTRODUCTION 

The binary system v O ctantis, consisting of a Kl III star 
ijHouk fc Cowlevl Il975l ) and a faint s econdary compone nt, 
identified as spectral type K7-M1 V l|Ramm et al.l [20081 ) is 
located in the southernmost portion of the southern hemi- 
sphere, i.e., i n close proxi mity to the Southern Celestial 
Pole. In 2009, iRamm et al.l found new detailed astrometric- 
spectroscopic orbital solutions, implying the possible exis- 
tence of a planelQ around the main stellar component. Thus, 
v Octantis A is one of the very few cases where a planet is 
found to orbit a giant rather than a late-type main-sequence 
star; the first case of this kind is the planetary-mass com- 
panion to the G9 III star HD 104985 (|Sato et alj|2003h . 

The composition of the v Octantis system leads to many 
open questions particularly about the orbital stability of the 
proposed planet. The semimajor axis of the binary system 
has been deduced as 2.55 ± 0.13 AU with an eccentricity 
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1 An alternative explanation involves a hierarchical triple system, 
invoking a p recessional motion of the main stellar component, as 
proposed bv lMorais fc Correial J2012I) . 



given as 0.2358 ± 0.0003 |Ramm et al.ll2009l ). Furthermore, 
the semimajor axis of the planet, with an estimated mini- 
mum mass MpSini of approximately 2.5 Mj, was derived as 
1.2 ± 0.1 AU (see Table [TJ. The distance parameters indi- 
cate that the proposed planet is expected be located about 
halfway between the two stellar components. This is a set- 
up in stark contrast to o ther observed binary systems host- 
ing one or more plane ts (jEggenberger. Udrv fc Mavoill2004l ; 
Ragh avan et al.ll200rj ) and, moreover, difficult to reconcile 
with standard orbital stability scenarios. 

A possible solution was proposed by lEberle fc Cuntzl 
|2010bl ). They argued that the planet might most likely be 
stable if assumed to be in a retrograde orbit relative to the 
motion of the secondary stellar component. Previous cases 
of planetary retrograde motion have been discovered, but 
with respect to the spin direction of the host stars. Exam- 
ples include (most likely) HAT-P-7b (or Kepler- 2b), WASP- 
8b an d WASP- 17b t hat w ere id entified by IWirm et all 
(2009), lOueloz et all |20ld ), and lAnderson et al. I |20ld h 
respectively. These discoveries provide significant chal- 
lenges to the standard paradigm of planet fo r matio n 
in protoplanetary disks as discussed by iLissaueJ (|l995t) . 
iLin. B odcnhci mer fc Richardsonl (1 19961 ) and others. 

The previous theoretical work by lEberle fc Cuntzl 
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<|2010bh explored two cases of orbital stability behaviour for 
v Octantis based on a fixed mass ratio for the stellar com- 
ponents and a starting distance ratio of pa = 0.379 for the 
putative planet (see definition below). They found that if 
the planet is placed in a prograde orbit about the stellar 
primary, it almost immediately faced orbital instability, as 
expected. However, if placed in a retrograde orbit, stability 
is encountered for at least 1 x 10 7 yr or 3.4 x 10 6 binary or- 
bits, i.e., the total time of simulation. This outcom e is also 
consis tent with the results from an earlier study by I Jeffervsl 
( 1974) based on a Henon stability analysis, which indicates 
significantly enlarged regions of stability for planets in ret- 
rograde orbits. 

Nonetheless, there i s a significant n e ed for augment- 
ing the earlier study by lEberle fc Cuntd (|2010bl ). Notable 
reasons include: first, lEberle fc Cuntz only considered one 
possible mass ratio and one possible set of values for the 
semimajor axis and eccentricity for the stellar components, 
thus disregarding their inherent observational uncertain- 
ties. Hence, it is unclear if or how the prograde and ret- 
rograde stability regions for the planetary motions will be 
affected if other a cceptable choice s of system parameters 
are made. Second, lEberle fc Cuntd chose the 9 o'clock po- 
sition as the starting position for the planetary component. 
Thus, it is unclear if or how the timescale of orbital sta- 
bility will change due to other positional choices in the 
view of previous studies, which have demonstrated a sig- 
nificant sensitivity in the outcome of stabili ty simulations 
on the adopted planetary starting an gle (e.g.. iFatuzzo et"al] 
20061; lYeager. Eberle fc Cuntd 1201 ll ). Third, and foremost, 



Eberle fc Cuntd did not utilize detailed stability criteria 



for the identification of the onset of orbital instability. 
Possible mathematical methods for treating the latter in- 
clude monitori ng the Jacobi constant and the zero velocity 
funct ion (e.g., Szcbchclvl 1 19671 ; lEberle. Cuntz fc Musielakl 
2008) as well as using the maximum Lyapun o v ex- 
ponent (e.g., iHilbornl 11994 ISmith fc Szebehelvl 1 19931 ; 
iLissauerl 1 19991 ; iMurrav fc Holmanl l200ll). In ou r prev i- 
ous work reported by ICuntz. Eberle fc Musielakl (2007). 
lEberle. Cuntz fc Musielakl (120081 ) and lQuarles et alj (|201ll i. 
we extensively used these methods for investigating the or- 
bital stability in the circular restricted 3-body problem. The 
same methods will be adopted in the present study. 

Our paper is structured as follows. In Sect. 2, we discuss 
our theoretical methods, including the initial conditions of 
the star-planet system and the stability criteria. The results 
and discussion, including detailed case studies, are given in 
Sect. 3. Finally, Sect. 4 conveys our conclusions. 



2 THEORETICAL APPROACH 
2.1 Basic Equations 

Using the orbital parameters attained by iRamm et al.l 
(2009), see Table [1] we consider the coplanar restricted 3- 
body problem (RTP) applied to the system of v Octantis. 
Following the standard conventions pertaining to the copla- 
nar RTP, we write the equations of motion in terms of the 
parameters p and po- The parameter p and the complemen- 
tary parameter a are defined by the two stellar masses mi 
and m2 (see below). The planetary distance ratio po de- 
pends on i?o and D, which denote the initial distance of 



the planet from its host star, the more massive of the two 
stars with mass mi, and the initial distance between the 
two stars, respectively. In addition, we represent the system 
in a rotating reference frame, which introduces Coriolis and 
centrifugal forces. The following eq uations of motion utilize 
these parameters l|Szebehelvll 19671 ): 
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The variables in the above equations describe how the 
state of a planet, commensurate to a test object, evolves due 
to the forces present. The state is represented in a Carte- 
sian coordinate system {x, y, z, u, v, w}. We use dot no- 
tation to represent the time derivatives of the coordinates 
|i = }• Since the variables {it, v, w} represent the ve- 
locities, the time derivatives of these variables represent the 
accelerations. The variables {n, r 2 } denote the distances of 
the planet from each star in the rotating reference frame. 
The value of / describes the true anomaly of the mass mi 
and e refers to the eccentricity of the stellar orbits. 



2.2 Initial Conditions 

The coplanar RTP indicates that the putative planet moves 
in the same orbital plane as the two stellar components; 
additionally, its mass is considered neglig iblfi The initial 
orbital velocity of the planet is calculated assuming a Ke- 
plerian circular orbit about its host star; see Fig. [1] for a 
depiction of the initial system configuration. For simplicity 
we assume the (u, w) components are equal to zero and the 
magnitude of the velocity to be in the direction of the v com- 
ponent. The mean value of the eccentricity of the binary, et, 
is used due to accuracy of the value given. A Taylor expan- 
sion in initial orbital velocity of the binary shows that the 
first order correction would only change the initial orbital ve- 
locity of the respective stars by a factor of 10~ 8 . The initial 
eccentricity of the planet, e p , has been chosen to be equal to 

2 Negligible mass means that although the body's motion is in- 
fluenced by the gravity of the two massive primaries, its mass is 
too low to affect the motions of the primaries. 



Stability in the v Octantis system 3 



Table 1. Stellar and Planetary Parameters of v Octantis 



Parameter 



Value 



Reference 



Spectral type (1) 
Spectral type (2) 
R.A. 
Decl. 

Distance (pc) 

My (1) 
My (2) 

Mi (M a ) 

M 2 (Mg) 
T eSA (K) 
Ri (Re) 
P b (d) 
a b (AU) 

Mpsim (Mj) 
a p (AU) 



Kl III 
K7-M1 V 
21 h 41 m 28.6463 s 
-77° 23' 24.167" 
21.20 ± 0.87 
2.10 ± 0.13 
~ 9.9 

1.4 ± 0.3 
0.5 ± 0.1 

4790 ± 105 
5.9 ± 0.4 
1050.11 ± 0.13 

2.55 ± 0.13 
0.2358 ± 0.0003 

2.5 ± ... 
1.2 ± 0.1 

0.123 ± 0.037 



Houk Cowley (1975) 



Ramm et al. (2009) 
ESA (1997) °''' 
ESA (l$9r) a ' b 
van Lceuwcn (2007) 
Ramm et al. (2009l c 
Drilling fc Landolt (20001 
Ramm et al. (20091 
Ramm et al. (20091 
Allende Prieto fc Lambert (19991 
Allende Prieto fc Lambert (19991 



Rammetal. 

Ramm et al. 

Rammetal. 
Ramrr^e^il^ 
Rammetal. 
Ramm et al. 



(2009) 
(20091 
(20091 
(2009) d 
(2009) d 
(2009) d 



" Data from SIMBAD, see [http://simabd.u-s trasbg.fr 
b Adopted from the Hipparcos catalog. 

c Derived from the stellar parallax, see lRamm et all ||2009| 1. 
d Controversial. 



zero due to the redundancy it provides in the selection of the 
initial planetary position. Noting that the inclination angle 
of the planet, i p , has not yet been determined, it is taken 
to be 90° consistent with the requirement of the coplanar 
RTP. 

To integrate the equations of mo tion, a sixth or der sym- 
plectic integration method is used IjYoshidal Il990l ). In the 
computation of Figs. [2] and [3] a fixed time step of 10 -3 is as- 
sumed. This proved to be appropriate since the average rela- 
tive error in the Jacobi constant is found to be on the order 
of 10 -8 . To increase the precision for the individual cases 
depicted in Figs. [S] to [5] a fixed time step of 10 -4 is used, 
which decreased the average relative error to the order of 



10' 



However, this choice entailed a considerable increase 



of the computation time. Therefore, the models pertaining 
to Fig. [2] were integrated for 10,000 binary orbits (approx- 
imately 29,000 years) and the models pertaining to Fig. [3] 
were integrated for 1,000 binary orbits (approximately 2,900 
years), noting that this figure is considerably less than uti- 
lized in the mainframe of our study for the determination of 
the statistical probability of a prograde planetary configu- 
ration. Moreover, the models of Figs. [5] to [9] were integrated 
for 80,000 binary orbits (approximately 232,000 years). 

In the treatment of this problem, four separate types of 
initial conditions of the system are considered. They encom- 
pass the two different starting positions of the planet with 
respect to the massive host star, which are: the 3 o'clock and 
9 o'clock positions (see Fig. [TJ. With each of these starting 
positions, we considered the possibility of planetary motion 
in counter-clockwise and clockwise direction relative to the 
direction of motion of the secondary stellar component; the 
directions are referred to as prograde and retrograde, re- 
spectively. The primary masses are considered to move in 
ellipses about the centre of mass of the system. 



2.3 Statistical Parameter Space 

In the present work we will explore the role(s) of small varia- 
tions in the system parameters. Notably we will investigate 
how the stellar mass ratio as well as the semimajor axis 
and inherently the eccentricity of the stellar components 
can change the outcome of our orbital stability simulations. 
Particularly when we consider the statistical probability of 
prograde orbits for the suggested planet. These variations 
are due to the observational uncer tainties of the respective 
parameter; see IRamm et all (|2009l ) and Tabled] for details. 

Equation (7) gives the definition of the mass ratio n 
with mi = 1.4±0.3 M Q and m 2 = 0.5±0.1 Mq. Taking mi 
and m2 as statistically independent from one another, the 
bounding values for fj, are [0.220, 0.306] , i.e., 0.263±0.043. In 
this case, we assumed the reported uncertainty bars for the 
stellar masses to be 1.25 a, noting that through this choice 
[i is consisten t with the mass ratio q — 0.38 ±0.03 derived by 
IRamm et all (2009) that was used for computing rai based 
on the sum of the two masses. Formulas and methods for the 
propagation-of-error, which in essence e mploy a first-o rder 
Tay lor expansion, have been given by, e.g.. I Me veil (1 19751 ); see 
alsoU ress et all l|l986l ) for numerical methods if the error is 
calculated via a Monte-Carlo simulation. 

Equation (12) gives the definition for the relative ini- 
tial distance of the planet, expressed in terms of Ro and D, 
i.e., the initial distance of the planet from its host star, the 
more massive of the two stars with mass mi and the initial 
distance between the two stars, respectively. This definition 
also allows to compute the b ounding values for pp. Fo llow- 
ing the assumption made by [Eberle fc Cuntj |2010bl ). the 
set up of the binary system will be initialized by assum- 
ing the secondary stellar component to be located in the 
apastron starting position, entailing that D is modified as 
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D = at (1 + £(,). Thus, the observation uncertainties of a p , 
Ob, and eb will be utilized to compute the bounding values 
for po (see Table [TJ. As result we obtain [0.344,0.418], i.e., 
0.381 ± 0.037. These values will be considered in the sub- 
sequent set of detailed orbital stability studies by making 
appropriate choices for fi and po when setting the initial 
conditions. The outcome of the various simulations will, in 
particular, allow a detailed statistical analysis of the mathe- 
matical probability for the existence of the suggested planet 
in a prograde orbit; see Sect. 13.31 

2.4 Stability Criteria 

In our approach we employ two independent criteria. One 
method deals with identifying events that can cause insta- 
bility, whereas the other method allows us to determine sta- 
bility. First, we monitor the relative error in the Jacobi con- 
stant to detect instabilities. By monitoring this constant 
we can determine when the integrator loses its accuracy 
or fails completely. It has been shown that this method is 
particularly sensitive to close approach events of the planet 
with any of the stell ar components. However, it has p revi- 
ously been shown by Eberle. Cuntz fc Musielakl |2008l ') and 
lEberle fe Cunt3 |2010aF that these close approach events are 
preceded by an encounter with the so-called "zero velocity 
contour" (ZVC) for the coplanar circular RTP. 

Generalizing this criteria to the elliptic RTP, an anal- 
ogous "zero velocity function" (ZVF) can be found that is, 
however, depend ent on the initia l distance between the stel- 
lar components ijSzebehelvll 19671). The ZVF can also be de- 
scribed by a pulsating ZVC; see ISzenkovits fc Makd HH) 
for updated results. Note that the ZVC will be largest when 
the binary components are at apastron and smallest when 
they are at periastron. Noting that we choose the stars to be 
at the apastron starting position, the largest ZVC is inher- 
ently defined as well. When the planetary body encounters 
the ZVC, the orbital velocity decreases dramatically, which 
reduces the u and v-related forces. As a result, the planetary 
body follows a trajectory that results in a close approach. 
This behaviour provides an indication that the system will 
eventually become unstable, resulting in an ejection of the 
planetary body from the system (most likely outcome) or 
its collision with one of the stellar components (less likely 
outcome) . 

The second c riterion utilizes the method of Lya- 
punov exponents l|Lvapunovl Il907l ). which is commonly 
used in nonlinear dynamics for determinin g the onset of 
chaos in different dynamic al systems (e.g.. iHilbornl fl994l ; 
iMusielak fc Musielakl l2009h . A numerical algorithm that is 
typically adopted for calculating the s pectrum of Lyap unov 
exponents was originally developed bv lWolf et al.l (| 1985h . In 
this approach, a set of tangent vectors Xi and their deriva- 
tives Xi with i = 1, . . . , 6 are introduced and initialized. In 
our approach, we choose all tangent vectors to be unit vec- 
tors for simplicity, and use a standard integrator akin to the 
Runge-Kutta schemes for computing changes in the tangent 
vectors within each timestep. After the tangent vectors be- 
come oriented along the flow, it is necessary to perform a 
Gram-Schmidt Renormalization (GSR) to orthogonalize the 
tangent space. In this specific procedure, we calculate a new 
orthonormal set of tangent vectors here denoted by primes 
(/), and obtain 



Xj 

x 2 



x 2 = 



X 2 , Xj 



X2 



X6 



(x 2 , X^x'J 



X<; 



X 6 , X 5 > X 5 



X 6 , Xi > Xi 



|x 6 - (x 6 , Xg)xg - (x 6 , x'^xij 



(13) 
(14) 

(15) 



With this new set of tangent vectors, the spectrum of 
Lyapunov exponents is computed by using 



A, 



A 4 _i + log 



3=1 



3-1 



(16) 



where A D = x D = 0. 

The Lyapunov exponents have extensively been used 
in studies of orbital stabili ty in different settings pertaining 
to the Solar Sys t em (e . g., Lecar, F ranklin fc Murison 1992; 
Milani fc Nobilil Il992l ; iLissaueii 1 19991 ; iMurrav fc Holmanl 
20011 ). Typically, the maximum Lyapunov exponent (MLE) 
is adopted as a measure of the rate with which nearby trajec- 
tories diverge. The reason is that the MLE offers the unique 
measure of the largest divergence. We previously used the 
MLE criterion to provide a cutoff for which stability can be 
assured given the behaviour of the MLE in t ime. This cut- 
off va lue has previously been determined by iQuarles et al.l 
|201ll) ; it is given as log A max = —0.82 on a logarithmic scale 
of base 10. 

Recalling that the cutoff in MLE is inversely related 
to the Lyapunov time th, the latter provides an estimate 
of the predictive period for which stability is guaranteed. 
The lower (i.e., more negative) the logarithmic value of the 
MLE, the longer the time for which orbital stability will be 
ensured. Due to the finite time of integration, this estimate 
is only valid for the integration timescale considered. But 
the nature of how the MLE changes in time will provide the 
evidence to whether the MLE is asymptotically decreasing 
and thereby allowing to identify orbital stability. Hence this 
cutoff will be used for establishing orbital stability of the 
suggested planet in the v Octantis system. 



3 RESULTS AND DISCUSSION 
3.1 Model Simulations 

In the following we report the general behaviour of sta- 
bility for the set of orbital stability simulations pursued 
in our study. Our initial focus is on the MLE. In Figs. [2] 
and [3] the various values of the MLE as a function of fj, and 
po are denoted using a colour code. The range of colours 
shows that dark red represents a relatively high (less nega- 
tive) logarithmic MLE value, whereas dark blue represents 
a relatively low (more negative) logarithmic MLE value. 
The white background represents conditions, which failed 
to complete the respective allotted run time; see Sect. 12.41 
for a discussion of possible conditions for terminating a sim- 
ulation. 

Figure [2] shows the result of our simulations for 10,000 
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Table 2. Survival Counts of the Suggested Planet in v Octantis for 1,000 binary orbits. 



Configuration 


Unstable Count 


Stable Count 


Stable Percentage 


Prograde 


3 o'clock 


17544 


1776 


9.2 % 


Prograde 


9 o'clock 


19200 


120 


0.6 % 


Retrograde 


3 o'clock 


1040 


19184 


94.9 % 


Retrograde 


9 o'clock 





20224 


100 % 



binary orbits, i.e., about 29,000 years, in the prograde con- 
figuration for starting conditions given by the 3 o'clock or 
9 o'clock planetary position. The results were obtained us- 
ing values of p, between 0.220 and 0.306 in increments of 
0.001. Concerning the initial conditions in po, denoting the 
relative starting position of the planet, we use values be- 
tween 0.272 and 0.502 in increments of 0.001. We show the 
lcr region in po between 0.344 and 0.418 where v Oct is pro- 
posed to exist by a grey shaded region; see Sect. 12.31 Between 
the two bounding values regarding po, there exists a region 
of stability as well as a marginal stability limit. The limit 
of marginal stability is con sistent with previous findings by 
iHolman fc Wiegert I (I1999T ) for the case with the planetary 
eccentricity taken into account. This region of marginal sta- 
bility is more pronounced in the 3 o'clock starting configu- 
ration; it begins at po ~ 0.315. The marginal stability region 
decreases as pi increases, which is expected because the per- 
turbing mass is becoming more significant as the mass ratio 
p, increases. There is a noticeable favoritism towards the 3 
o'clock starting position as indicated by the survival counts 
given in Table[2] This outcome is expected as the perturbing 
mass exerts a lesser force on the planet if initially placed at 
the 3 o'clock position. 

Figure [3] shows the result for the retrograde configura- 
tion for 1,000 binary orbits considering both the 3 o'clock 
and 9 o'clock planetary starting positions. Due to the larger 
number of conditions that survived during this time interval 
we have reduced the bo unds of initial c o ndition s following 
the bounding values of lEberle fc Cunt3 (|2010bl ). Here the 
bounding values for p are given as 0.2593 and 0.2908 with 
increments of 0.0001. The increment in pi was reduced to 
allow better insight into the role of chaos in this type of 
configuration. The initial conditions of po begin with 0.353 
and end with 0.416 in increments of 0.001, which allows to 
cover the most pronounced features of the retrograde con- 
figuration. 

Figure [3] demonstrates a broad region of stability within 
the selected parameter space. The 9 o'clock starting posi- 
tion has many regions of bounded chaos. These regions of 
bounded chaos exhibit locally less negative values for MLE, 
which however do not exceed the previously determined cut- 
off value of log A max = —0.82; see Sect. 12.41 The nature of 
these orbits and the connection to the MLE will be dis- 
cussed in Sect. 13.21 The greatest region of bounded chaos 
is obtained for po sa 0.41. This region decreases in po as p 
increases for the same reason as the change in the marginal 
stability limit in Fig. [2] The 3 o'clock starting position shows 
similar features but there exists a region of instability with 
this starting position at po ~ 0.41. This region coincides 
with the extended region of bounded chaos in the 9 o'clock 



starting position. The nature of this region c an be deduced 
from the parameters given bv lRamm et al.l (|2009l ); see Ta- 
ble [1] There exists a 5:2 resonance of the planet with the 
binary that places the suggested planet moving towards the 
barycentre at the periastron point of the binary inducing an 
instability. For the 9 o'clock position the suggested planet 
would be moving away which would provide a large pertur- 
bation but insufficient for producing an instability. 

Figure [4] provides the histograms of simulations per- 
taining to both starting positions for the prograde configu- 
ration. These histograms indicate the number of surviving 
cases related to the array of initial conditions for the pa- 
rameter space detailed in Fig. [5] The bin width is 0.01 in 
the logarithmic MLE scale, corresponding to 1.47 x 10~ 4 in 
the linear MLE scale. Although there is not a case that sur- 
vives 10,000 binary orbits in the 1 a region, we provide these 
histograms to demonstrate what different initial conditions 
entail near the lower and upper uncertainty limits. Figure [2] 
shows that all the surviving cases cluster towards the lower 
bound in po. Furthermore, it is obtained that the 9 o'clock 
starting position has far fewer surviving initial conditions 
than the 3 o'clock starting position. This is also reflected 
in the histograms. The histogram of the 9 o'clock position 
is sparse; without the aid of the 3 o'clock histogram there 
would be too few counts to provide substantial information. 
But inspecting the 3 o'clock histogram reveals that there 
is a Gaussian-type distribution about the logarithmic MLE 
value of —2.2. This indicates that for the limited simulation 
time adopted for this part of our study, there is a very small 
number of cases indicating stability for the prograde con- 
figuration. However, these cases are near the 2cr or 3er limit 
concerning the permitted values of the planetary position po 
as discussed in Sect. 13.31 

3.2 Case Studies 

In the following we probe some orbits of retrograde configu- 
rations in more detail; see Figs.[5]to[9]for information. These 
figures show from left to right the respective orbital diagram 
in a rotating reference frame, the associated Lyapunov spec- 
trum, and the Fourier periodogram. Through these diagrams 
we can assess the nature of these orbits to determine stabil- 
ity and measure the chaos within the system. 

Figure [5] with (p, p ) = (0.2825, 0.400) is provided to 
show a case of instability. In this case we find that the orbit 
presents itself as a wide annulus and already at first glance 
appears chaotic. The Lyapunov spectrum reveals that the 
MLE levels off at log A = — 1 fairly quickly, which is near 
the proposed cutoff value for unstable chaotic behaviour. For 
the timespan between 300 and 800 binary orbits, the MLE 
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Table 3. Statistical Probability of Prograde Planetary Orbits". 



Mass Ratio 9 o'clock Position 3 o'clock Position 



p a Probability (%) a Probability (%) 



0.306 ... ... 2.30 1.1 

0.263 ... ... 2.11 1.7 

0.220 2.95 0.16 1.89 2.9 



Based on 10,000 binary orbits; see Fig. [2] 

fluctuates until it finally pushes upward to A max . The Fourier 
periodogram shows that this case is truly chaotic due to the 
amount of noise presented by the periodogram. There are 
many periods which may pr omote the outcome of i nstability 
via r esonance overlap (e.g.. iMudrvk fe Wu 2006:; iMardlind 
I2007T ). The proposed planet is approaching the barycentre 
when the binary is at its periastron position, which should 
be considered the main cause of this instability. 

Figures [6] and [7] with (p,,po) = (0.2805,0.358) and 
(0.2630, 0.374), respectively, exhibit similar features between 
each other in the orbit diagram and periodogram. However, 
there is a difference in the respective Lyapunov spectra. Both 
cases show trends of stability, although Fig. [6] demonstrates 
bounded chaos. It presents a different nature than the pre- 
vious case of instability (see Fig. [5| as the MLE settles near 
log A = —1.8 and does not fluctuate enough to reach A ma x- 
Both cases given as Figs. |S] and [7] represent orbits that form 
annuli; however, the annulus widths are less than that of 
the unstable case depicted in Fig. [S] They are thus classi- 
fied as stable with Fig. [5] being chaotic and Fig. [7] being 
non-chaotic. 

Figures H and H with (p, p ) = (0.2696, 0.401) and 
(0.2780,0.388), respectively, reveal similar features com- 
pared to the previous set of figures regarding all three di- 
agrams. These cases are provided to demonstrate the simi- 
larities between the 3 o'clock starting positions depicted in 
Figs. [6] and [7] and the 9 o'clock starting positions depicted 
in Figs. [5] and [5] The comparison between Figs. |S] and [5] in- 
dicates a similar orbit diagram and periodogram, but the 
limiting values of the MLE settle at a slightly more negative 
value in Fig. HJ] Figures [7] and [5] supply similar orbital dia- 
grams but there are minor differences in the periodograms 
and Lypaunov spectra. These figures display the same be- 
haviour in MLE but there is a difference in the periodic 
nature of the decrease in MLE, which is also revealed by 
the difference in the number of peaks in the corresponding 
periodograms. 

3.3 Statistical Analysis of Orbital Stability 

A significant component of the present study is to provide a 
statistical analysis regarding whether or not the suggested 
planet in the v Octantis system is able to exist in a pro- 
grade or retrograde orbit. Previously, we reported on the 
general behaviour of orbital stability for the set of simula- 
tions with focus on the MLE (see Sect. 12. 4j) . In Figs. [2] and [3j 
the various values of the MLE were given as function of the 



mass ratio p and the initial planetary distance parameter po. 
The bounding values for p were given as [0.220, 0.306], i.e., 

0. 263 ±0.043, and for p as [0.344,0.418], i.e., 0.381 ±0.037. 
We found that based on the MLE as well as a number count 
of surviving cases (see Sect. 13. ip that the existence of the 
suggested planet in a retrograde orbit is almost certainly 
possible, whereas in a prograde orbit it is virtually impossi- 
ble. 

For prograde orbits based on the "landscape of survival" 
(see Fig. 0, another type of analysis can be given, see Ta- 
ble [3] We conclude that for the 9 o'clock planetary starting 
position and for a relatively high stellar mass ratios, i.e., 
p > 0.25, the planet must start very close to the star, i.e., 
that the associated probability for po is way outside the sta- 
tistical 3<t limit. Only for relatively small stellar mass ratios, 

1. e., p ~ 0.22, the planetary starting distance po would be 
barely consistent with the 3a limit, amounting to a statisti- 
cal probability of 0.16%. For a 3 o'clock planetary starting 
position, these probabilities pertaining to the permitted val- 
ues of po are moderately increased. It is found that the as- 
sociated probabilities for p given as 0.306, 0.263, and 0.220 
are 1.1%, 1.7%, and 2.9%, respectively. 

Moreover, it should be noted that this statistical anal- 
ysis refers to a "landscape of survival" that has been ob- 
tained for a timespan of 10,000 binary orbits, corresponding 
to 29,000 years. An analysis of the evolution of that land- 
scape shows that it tends to dissipate if larger timespans are 
adopted. Even for timespans of 10,000 binary orbits, there 
are many parameter combinations (p, po) where, concerning 
prograde orbits, no orbital stability is found even if p values 
are chosen at or beyond the statistical 3<r level. 

Another perspective is offered through Fig. 1101 which 
provides additional insight into the statistical nature of our 
simulations. It addresses the view of long term stability. Fig- 
ure [10] depicts two curves representing the number of sur- 
viving configurations as a function of time for each starting 
position in the prograde configuration. The features of this 
logarithmic plot show us that for the first 10 binary orbits 
there is not a preference towards stability of either starting 
position. After 10 binary orbits the 9 o'clock position shows 
a significantly steeper trend of instability as indicated by 
the slope of the respective curve. The 3 o'clock position also 
shows a trend toward instability but it is asymptotically ap- 
proaching a limit of 10 3 much more smoothly. With these 
limits, we can estimate the likelihood for a prograde plane- 
tary configuration to be extremely small in the framework 
of long-term simulations. However, for a very large range of 
(p, po) combinations, the stability of retrograde planetary 
configuration appears almost certainly guaranteed, even for 
timescales beyond 10 5 binary orbits. 



4 CONCLUSIONS 

This study offers detailed insights into the principal possi- 
bility of prograde and retrograde orbits for the suggested 
planet in the v Octantis system, if confirmed through future 
observations. A previous study by lEberle fc Cuntzl i|2010bh 
concluded that the planet is most likely stable if assumed to 
be in a retrograde orbit with respect to the secondary sys- 
tem component. In the present work, we were able to confirm 
this theoretical finding, while taking into account the obser- 
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vationally deduced uncertainty ranges of the orbital param- 
eters for the stellar components and the suggested planet as 
well as different mass ratios of the stellar components. 

In addition, we also employed additional mathematical 
methods, particularly the behaviour of the maximum Lya- 
punov exponent (MLE). We found that virtually all cases 
of prograde orbit became unstable over a short amount of 
time; the best cases for survival (although barely significant 
in a statistical framework) were found if a relatively small 
mass ratio for the two stellar components was assumed and 
if the initial configuration was set for the planet to start at 
the 3 o'clock position. This later preference is also consistent 
with previous general results bv lHolman & Wiegert I (|l999l ) 
and others. 

Thus, our results of the retrograde simulations indicate 
a high probability for stability within this type of configu- 
ration. Moreover, this configuration displays regions in the 
system where bounded chaos can exist due to the resonances 
present. Concerning retrograde orbits, there is little prefer- 
ence for the planetary starting position. In summary, we 
conclude that the v Octantis star-planet system is an inter- 
esting case for future observational and theoretical studies, 
including additional long-term orbital stability analyses if 
the suggested planet is confirmed by follow-up observations. 
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X* (AU) 

Figure 1. Example of the possible initial starting configurations. The directions of the initial velocities of the stellar components (y Oct A 
and v Oct B) are shown with black coloured arrows and appropriate labels. The proposed planet is represented by the red coloured dots 
in the 3 o'clock position (right) and the 9 o'clock position (left). The possible directions of the initial velocity for the planetary component 
are either prograde (purple) or retrograde (blue). The arrows indicate the directions in the initial configurations of the system. 
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Po Po 

Figure 2. Detailed results of simulations for the parameter space for the case of prograde planetary motion. The left and right panel 
represent the simulations with the planet initially placed at the 9 o'clock and 3 o'clock position, respectively. The grey area represents 
the statistical la range in po, which is 0.381 ± 0.037. 
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Full Space Histogram of Nu Oct at 9 o clock starting position 
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Figure 4. Histograms of configurations that survived 1000 binary orbits in prograde motion. Note the significant difference in scale 
regarding both the x and y-axis. 
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Figure 5. Case study of planetary motion with the planet placed in the 3 o'clock position and in retrograde motion. This case displays 
the conditions for fi = 0.2825 and po = 0.400. 
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Figure 6. Case study of planetary motion with the planet placed in the 3 o'clock position and in retrograde motion. This case displays 
the conditions for fi = 0.2805 and po = 0.358. 
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Figure 7. Case study of planetary motion with the planet placed in the 3 o'clock position and in retrograde motion. This case displays 
the conditions for fi = 0.2630 and po = 0.374. 
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Figure 8. Case study of planetary motion with the planet placed in the 9 o'clock position and in retrograde motion. This case displays 
the conditions for fi = 0.2696 and po = 0.401. 
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Figure 9. Case study of planetary motion with the planet placed in the 9 o'clock position and in retrograde motion. This case displays 
the conditions for = 0.2780 and p = 0.388. 
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Figure 10. Logarithmic representation of the surviving number of configurations (prograde orbits) as function of time given in units of 
completed binary orbits. 



